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One contribution to the energy loss of hard partons propagating through a medium as created in 
ultrarelativistic heavy-ion (A-A) collisions are elastic scattering processes with medium constituents. 
The magnitude of this energy loss contribution depends crucially on the effective mass of the medium 
constituents — in the limit of non-recoiling static scattering centers, elastic energy loss vanishes. 
Thus, it is important to constrain the amount of elastic energy loss in order to gain information 
about the nature of the degrees of freedom in the medium as resolved by a hard parton. So far, 
the relative fraction of elastic (or rather incoherent) energy loss has been constrained from above 
by using pathlength dependent observables. However, using the observation that subleading gluon 
energy dissipation into the medium is probed by some observables, such as hard dihadron correlations 
or the dijet asymmetry, a constraint from below can also be found. In this work, this idea is worked 
out within the in-medium shower evolution Monte Carlo (MC) code YaJEM-D. 

PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 

The energy loss of hard partons propagating through 
the soft medium created in heavy-ion collisions, leading 
to a suppression of high Pt hadron spectra, has long 
been regarded as a promising tool to gain information on 
medium properties [ll-Q- However, there is a persistent 
puzzle in the literature with regard to the balance be- 
tween radiative energy loss (where energy is carried by 
medium-induced gluon radiation) and elastic energy loss 
(where energy is carried in the recoil of medium scatter- 
ing centers). 

If one starts with a first principles calculation from 
perturbative Quantum Chromodynamics (pQCD), mod- 
elling the medium as thermal distributions of quasi-free 
quarks and gluons with thermal masses of order gT 
(where g is the strong coupling with as — (7^/(47r) and T 
the medium temperature) as motivated by thermal field 
theory, one naturally arrives at the conclusion that elastic 
energy loss is large |7H14| , possibly capable of accounting 
for about half of the energy carried away from a light 
quark passing through a medium. 

If, on the other hand, one starts with the observa- 
tion that elastic energy loss is an incoherent process and 
hence has a linear pathlength dependence in a constant 
medium, one can analyze the experimentally measured 
pathlength dependence by comparing models with the 
nuclear suppression factor Raa in non-central collisions 
as a function of the angle (j) with the reaction plane. Such 
comparisons, both for a generic phenomenological model 
[Tsj and a detailed Monte-Carlo (MC) model [H] find 
that incoherent contributions to the total energy loss 
must be small and of the order of about 10%. To the 
degree that elastic energy loss is only a part of all inco- 
herent processes, this severely constrains the amount of 
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elastic energy loss taking place in nature. 

The inevitable conclusion is that the main assumption 
made in pQCD calculations of elastic energy loss, i.e. 
that the medium DOF are almost free (quasi)-particles 
which can take a sizeable amount of recoil energy away 
from a leading parton does not appear to be true in na- 
ture. In this case, constraining the amount of elastic 
energy loss offers information as to the nature of the scat- 
tering centers in the medium. 

In this paper, we investigate a way to constrain elas- 
tic energy loss different from pathlength dependence by 
studying the interaction of subleading jet fragments. The 
mechanism is the same which has been invoked to ex- 
plain the dijet asymmetry observed by ATLAS and CMS 
[ITI fisj — low Pt shower partons are likely to be scat- 
tered out of a jet cone, and hence the combination of ra- 
diative energy loss as a source of soft gluons and elastic 
interactions decorrelating those from the jet can plausi- 
bly account for the observed dijet asymmetry jioj. Here, 
we investigate the same effect in the context of triggered 
back-to-back correlations, which offer a somewhat less 
complicated environment free from the problems of jet 
reconstruction in a medium [201 . 



II. THE MODEL 

A. General framework 

We compute the strength of back-to-back hadron 
correlations by combining a MCpQCD calculation of 
back-to-back parton production |2ll . |22| with (geometry- 
dependent) medium-modified fragmentation functions 
obtained with the MC code YaJEM |2i,|2l] with the min- 
imum virtuality scale down to which the shower is evolved 
in the medium determined by the prescription outlined in 
[25| (YaJEM-D). Note that our modelling is constrained 
by multiple observables, among them the reaction plane 
angle dependence of the nuclear suppression factor Raa 
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for non-central collisions [25| at RHIC, the reaction plane 
dependence of the dihadron correlation suppression fac- 
tor Iaa [13 at RHIC and the nuclear suppression at LHC 
[26j , thus we are only allowed to modify our modelling 
in a way that does not destroy agreement with these ob- 
servables. 

In LO pQCD, the production of two hard partons fc, I 
is described by 



da 



AB^kl+X 



dpj.dyidy2 



= >Ja;i/,/^(a;i,Q^)a;2/j/B(a;2,Q^) 



dt 



(1) 

where A and B stand for the colliding objects (protons 
or nuclei) and 2/1(2) is the rapidity of parton k{l). The 
distribution function of a parton type i in A at a mo- 
mentum fraction xi and a factorization scale Q ^ px is 
fi/A{xi,Q^). The distribution functions are different for 
free protons [13, [13 and nucleons in nuclei [H-[3lj . The 
fractional momenta of the colliding partons i, j are given 
by xi^2 = ^ (exp[±2/i] -|- exp[±y2])- Expressions for the 

pQCD subprocesses ^- — {s,i,u) as a function of the 

parton Mandelstam variables s, t and it can be found e.g. 
in [32]. 

To account for various effects, including higher order 
pQCD radiation, transverse motion of partons in the nu- 
cleon (nuclear) wave function and effectively also the fact 
that hadronization is not a collinear process, the distri- 
bution is commonly folded with an intrinsic transverse 
momentum kx with a Gaussian distribution, thus creat- 
ing a momentum imbalance between the two partons as 

PTi + PT2 = kT- 

Eq. ([T]) is evaluated at midrapidity yi — y2 — and 
sampled using a MC code introduced in [2l| by first gen- 
erating the momentum scale of the pair and then the 
(momentum-dependent) identity of the partons. A ran- 
domly chosen kx with a Gaussian distribution of width 
2.5 GeV is then added to the pair momentum. 

Under the assumption that the distribution of vertices 
follows binary collision scaling as appropriate for a LO 
pQCD calculation, the probability density to find a ver- 
tex in the transverse plane is 



Pixo,yo) = 



TA{ro + h/2)TA{ro~h/2) 

TAA{h) 



(2) 



where the thickness function is given in terms of Woods- 
Saxon distributions of the the nuclear density /9^(r, z) 
as ^^(r) = J dzpA{r,z) and TaaO^) is the standard nu- 
clear overlap function Taa{^) = J c^^s Ta(s)Ta(s — b) 
for impact parameter b. Each parton pair is placed at a 
probabilistically chosen vertex (xq, yo) sampled from this 
distribution with a random orientation (jj with respect to 
the reaction plane. 

Both partons are then propagated on eikonal paths 
through a hydrodynamical medium ^SJj and for this path 
the leading and first two subleading fragments are com- 
puted using the medium- modified conditional probability 



densities Ai{zi,fi), ^2(^1, 2^2, A^) and A3{zi, Z2, z^, fi) sa 
A2{zi + Z2, Z3, ^) given the path [33| . For this, we utilize 
YaJEM-D. Finally we check if there is a hadron in the 
event which fulfills the trigger condition. If not, we dis- 
card the event and start generating a new one. If there 
is a trigger hadron, we bin the remaining hadrons in the 
event on the near and away side in either p^^'^°'^ or zt- 



B. Medium-modified fragmentation 

The key ingredient containing the information about 
the medium evolution and shower-medium interaction is 
the medium-modified fragmentation function (MMFF) 
which wc compute in YaJEM-D. The MC code YaJEM- 
D is based on the PYSHOW code [35^ which is part of 
PYTHIA [36,]. It simulates the evolution from a highly 
virtual initial parton to a shower of partons at lower vir- 
tuality in the presence of a medium down to a minimum 
scale Qo ~ \/E/L where E is the energy of the shower 
initiator and L is the in-medium pathlength. A detailed 
description of the model can be found in |23l - l25j ] . 

In the RAD (radiative energy loss) scenario, a parton 
a gains virtuality during its propagation time Tq as 



(3) 



which leads to extra branchings and soft gluon pro- 
duction where g(C) is the local rate of virtuality gain. 
However, YaJEM-D also implements the so-called DRAG 
scenario in which propagating partons lose energy (and 
momentum) according to 



dCDiC) 



(4) 



where -D(C) is the local drag coefficient correspond- 
ing to a mean energy loss dE /dC, per unit length. This 
model is also suited to treat incoherent, elastic energy 
loss. These coefficients are tied to the hydrodynamical 
parameters of the medium evolution model via 



q[D]{0 = K[Kd] ■ 2 • [e(C)]3/^(coshp(C) - sinhp(C) cos^) 

(5) 

where V' is the angle between bulk medium flow and 
the parton direction, p is the flow rapidity and e is the 
medium energy density. This expression is a proxy for 
the Lorentz-contracted density of scattering centers in 
the fiowing medium as seen from a moving parton. K and 
Kjj then parametrize the actual strength of the parton- 
medium interaction given this density. 

In |2^] it was found that both scenarios obey a scal- 
ing law, i.e. to a good approximation what matters for 
the end result are the integrals /S.E = J d(D{C,) and 
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AQ^ = J d(q{() rather than the fuh functional depen- 
dence of the coefficients on C. It was also found that 
both scenarios give equally good descriptions of the nu- 
clear suppression factor Raa at RHIC kinematics (their 
results do, however, differ for other more differential ob- 
servables such as jet shapes [131 )• This implies that suit- 
able combinations of ((j, D) where one of the parameters 
does not vanish will also describe the Pt dependence of 
Raa- 
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FIG. 1: (Color online) A comparison of medium-modified 
fragmentation functions for a 30 GeV quark computed in 
YaJEM-D for various relative contributions of radiative and 
elastic energy loss, indicated by the mean accumulated vir- 
tuality AQ^ or the mean energy loss AE along the parton 
path. 

We can gain some insight into this by studying the 
MMFF for different parameter combinations {q,D) (or 
rather their path-integrated values (AQ^, AE)) in Fig.[T] 
Folded with a pQCD parton spectrum, one finds that this 
fragmentation function has a mean z of ~ 0.3, i.e. just 
where all scenarios agree. This explains why there is 
no pronounced difference between the scenarios for sin- 
gle inclusive hadron production in spite of the obvious 
differences in the MMFF at high and low z. 

It is difficult to exploit the differences at high z to 
determine which behaviour of the MMFF is realized in 
nature since high z fragmentation is a rare process. How- 
ever, the low z behaviour is easily studied by observing 
subleading shower hadrons in back-to-back dihadron cor- 
relations on the away side. 



III. RESULTS 

We show the away side hard dihadron suppression fac- 
tor Iaa as a function of zt (momentum of observed away 
side hadron divided by trigger hadron momentum) com- 
puted as outlined above in Fig. [5] and compare with 
STAR data [11]. We do the comparison for two dif- 
ferent scenarios: One (YaJEM-D) assumes pure radia- 
tive energy loss (and agrees with the result presented in 



[111), the other (YaJEM-DE) uses q = 0.8 ■ qmax and 
D = 0.1 • Dmax (where (jmax and Dmax are the param- 
eter values one extracts for a fit to Raa assuming that 
the corresponding other parameter is zero, i.e. for a pure 
radiative or a pure elastic energy loss scenario). These 
values are motivated by the constraint that the elastic 
contribution should be no more than of order 10% of 
the total energy loss as constrained by pathlength de- 
pendence. They also satisfy the constraint that the Pt 
dependence of single hadron Raa should be unchanged 
within experimental errors. 
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FIG. 2: (Color online) The away side Iaa for RHIC Au-Au 
collisions, computed using the in-medium shower MC code 
YaJEM-D with only radiative energy loss or with a 10% con- 
tribution of elastic energy loss (YaJEM-DE) compared with 
STAR data [H. 

We see that at high zt there is no statistically sig- 
nificant difference between the two scenarios. However, 
at low Zt they clearly separate, in accordance with the 
behaviour seen in the MMFF in Fig. [1] and the scenario 
including elastic energy loss describes the data much bet- 
ter in this region. 

The physics interpretation of this result is simple: 
Medium-induced radiative energy loss creates a large 
number of soft gluons, which in turn leads not only to a 
depletion of the MMFF at high z ('quenching') but also 
to a rapid growth at low z and consequently of I'aa^ at 
low Zt- a certain contribution of elastic energy loss is 
needed to dissipate the energy of these soft gluons into 
the medium and thus tame the growth of the MMFF. 

While it is difficult to extract relative fractions of en- 
ergy loss in a medium-modified shower as such since it is 
impossible to tag any particular branching and assign it 
to either vacuum radiation or medium-induced radiation 
(see discussion in '23] ), based on the magnitude of the 
transport coefficients compared to pure radiative or elas- 
tic scenarios, it seems that a fraction of about 10% elastic 
contribution which is allowed by pathlength dependence 
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is also sufficient to account for the observed modification 
of the fragmentation function in dihadron correlations. 
However, note that while it is true that choosing a value 
of 10% D^ax corresponds to also 10% mean energy loss 
(the actual energy loss in any single MC event may dif- 
fer) with respect to the value obtained with D,nax^ the 
same relation is not true for q which has a complicated 
non-linear relation with gluon radiation which includes 
LPM suppression of soft gluons during decoherence time 
and kinematical restrictions. Thus, ratios of the coeffi- 
cients D,q do not imply the same ratio of mean energy 
loss from the leading parton. Therefore, a more precise 
whay of thinking about the problem is fixing the relative 
strength of the parameters {K^ Kb) which determine the 
value of (g, D) given the hydrodynamical parameters of 
the medium evolution model. 

Given that the same mechanism, i.e. dissipation of 
soft gluon energy into the medium, is used to explain 
the diiet asymmetry [li| observed by ATLAS and CMS 
[13, US], a crucial test for the combination of radiative and 
elastic energy loss determined in the scenario YaJEM- 
DE is whether or not the model is able to account for 



the measured dijet asymmetry. This question will be 
addressed in future work. 

If the idea that elastic energy loss is observed dissi- 
pating the energy of low z gluons turns out to be cor- 
rect, pathlength dependence can be used to constrain the 
amount of incoherent energy loss from above whereas dis- 
sipation of energy in soft gluons can be used to constrain 
it from below, i.e. a relatively precise determination of 
incoherent energy loss should be viable. This in turn may 
greatly help to address the question what the microscop- 
ical degrees of freedom in the medium as probed by a 
hard parton are. 
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